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ABSTRACT : Results are presented of parametric studies to assess the effects of various 
parameters on the free vibration behavior (natural frequencies) of [±0] 2 angle-ply fiber 
composite thin shells in a hot environment. These results were obtained by using a three- 
dimensional finite element structural analysis computer code. The fiber composite shell is 
assumed to be cylindrical and made from T300 graphite fibers embedded in an 
intermediate-modulus high-strength matrix (IMHS). The residual stresses induced into 
the laminated structure during the curing are taken into account. The following parameters 
are investigated: the length and the thickness of the shell, the fiber orientations, the fiber 
volume fraction, the temperature profile through the thickness of the laminate and the dif- 
ferent ply thicknesses. Results obtained indicate that: the fiber orientations and the length 
of the laminated shell had significant effect on the natural frequencies. The fiber volume 
fraction, the laminate thickness and the temperature profile through the shell thickness had 
a weak effect on the natural frequencies. Finally, the laminates with different ply thick- 
nesses had insignificant influence on the behavior of the vibrated laminated shell. 

KEY WORDS: laminated cylinders, fiber composites, angle-ply laminates, composite 
structures, computational simulation, residual stresses, high temperature, moisture, struc- 
tural analysis, finite element analysis, free vibration, natural frequencies, mode shape. 

INTRODUCTION 

H igh speed flight vehicles under hostile environments at speeds which are 
multiples of the speed of sound require new materials and improved compu- 
tational approaches to design them. New materials, known as elevated tempera- 
ture composites, have high strength, are light in weight, and can be tailored for 
specific performance requirements. These materials are already finding applica- 
tions and acceptance in aircraft frame and engine structures. Cost effective usage 
of elevated temperature composite materials requires new ideas and innovative 
concepts in the areas of material selection and structure fabrication in order to 
design for durability by properly accounting for many failure mechanisms in- 
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herent in composites. The influence of hostile environments, nonlinear material 
and structural behavior, and the coupling between responses induced by various 
loads require sophisticated analysis methods. 

A stand-alone multi-disciplinary computer code, CSTEM (Coupled Structural/ 
Thermal/Electro-Magnetic Analysis/Tailoring) [1 to 5] has been developed by in- 
tegrating three-dimensional finite element structural/stress analysis methods with 
several single-discipline codes including those for integrated composite 
mechanics [6,7]. 

This paper demonstrates the use of CSTEM to computationally simulate air- 
craft engine composite structures subjected to free vibration at elevated tempera- 
tures. Knowledge of the natural frequencies of structures is important in consid- 
ering their response to various kinds of excitation especially when structures and 
force systems are complex and when excitations are not periodic [8,9]. The struc- 
ture studied is a thin, cylindrical composite shell. The composite is to be made 
from angle-ply [ ± 0} 2 laminates containing plies with uniform thickness. The re- 
sidual stresses induced into the laminated composite during the curing process 
are taken into account. Parametric studies were performed to examine the effects 
on the natural frequencies of parameters such as: the length and the laminate 
thickness of the shell, the temperature profile through the thickness of the lami- 
nate, the fiber volume fraction, the fiber orientations and laminates with different 
ply thicknesses. 

BRIEF DESCRIPTION OF THE SIMULATION PROCEDURE 

The general purpose computer program CSTEM |1] computationally simulates 
the coupled multidisciplinary structural, heat transfer, vibration, acoustic, and 
electromagnetic behavior of elevated temperature, layered multi-material com- 
posite structures in aggressive environments. All the disciplines are coupled for 
nonlinear geometrical, material, loading, and environmental effects. CSTEM is 
based on a modular structure, (Figure 1), and is accessed by the user through its 
executive module. The structural response can be predicted at all the composite 
scales including constituents (fiber and matrix) and ply level. 

The composite mechanics is simulated with the ICAN (Integrated Composite 
Analysis) module [6,7], which performs, among others, through-the-thickness 
point stress analysis. The micro-mechanics equations embedded in ICAN in- 
clude the effects of temperature and moisture. ICAN simulates the behavior of 
polymer composites from the constituent material level to the laminate level as 
depicted in Figures 2 and 3. ICAN includes a resident material property data 
bank for commercially available fiber and matrix constituent materials at room 
temperature. The user needs only to specify a code name for the desired constitu- 
ent materials (rather then having to manually input all the properties) in the 
CSTEM input. 

The procedure to computationally simulate composite structures is shown in 
Figure 2 and consists of four parts: the constituents, the synthesis, the finite ele- 
ment structural analysis, and the decomposition. For detailed description of these 
four parts see References [6] and [7] . The computational procedure steps are: 
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Figure 1 . CSTEM modular structure. 


Compone nt 



Figure 2. Integrated composite analysis. 
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Figure 3. Regions of constituent materials and nonlinear material characterization model. 
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1. Constituents -The material properties of the matrix, in any operational tem- 
perature and moisture condition are updated based on the constitutive equa- 
tion depicted in Figure 3. 

2. Synthesis —The thermal and mechanical properties of the plies are synthesized 
knowing the mechanical properties of the fibers and the matrix using compos- 
ite micro-mechanics theory. 

The extensional stiffness matrix A, the coupling stiffness matrix B, and the 
bending stiffness matrix D are computed by using the laminate theory as well 
as the effective laminate thermal and mechanical engineering properties [10], 

3. Finite Element Analysis - The finite element analysis is performed at the com- 
posite structure scale, and among the others, the resulting forces and mo- 
ments are computed at each node of the finite element mesh. Because ICAN 
performs through -the -thickness point stress analysis, steps (1) and (2), as well 
as (4) are performed at each individual node of the finite element mesh. 

4. Decomposition— The stresses and strains in each ply are computed by using 
the laminate theory. The stresses and strains in the matrix and the fibers for 
each ply are computed by using composite micromechanics. 

The CSTEM code is written in the FORTRAN programming language and in- 
stalled on the NASA Lewis CRAYYMP computer system. 

SUMMARY OF THE BASIC EQUATIONS 
Ply Level Equations 

It is assumed that the heterogeneous ply consisting of fibrous and matrix 
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phases can be modeled as a homogeneous and orthotropic medium. Also, the ply 
is assumed to be subjected to a state of plane stress. The stress strain relationship 
for the kth ply in material axes 1 , 2, 3 (1 is the direction of the fibers and 3 is the 
direction normal to the 1,2 plane) is the following: 

l^/i = [0/.l(k/l — — k/dw) (0 


where 


\oiV = l ff,,a 2 ,a 12 ]. 


1 0i. I 


G,. 0,2 o 

021 022 0 

0 0 033 


and 0 u = EJP, Q 21 = E 2 /P , 0 12 = 0 21 = * 12 EJP = v 2x E 2 IP y 0 33 = G x2 
and P = 1 ~ v x2 v lx . The effective properties of the kth ply, E 2y v X i, v lx and 
G ,2 are computed analytically from the constituent material properties (the fiber 
and the matrix) [6,7], 



a, and a 2 = 

0 i and 0 2 = 

AT = 
AM = 


|f, £ 27 , 2 ] is the mechanical strain vector 
I a x a 2 0 ]AT is the thermal strain vector and 
[0, 0 2 O]AAf is the moisture strain vector 

the effective free thermal expansion coefficients of the kth ply with 
respect to the 1 - and 2 -axes 

the effect free moisture expansion coefficients of the Ath ply with 

respect to 1 - and 2 -axes 

the temperature in the kth ply 

the moisture in the Arth ply. 


From Equation ( 1 ) we have: 


\Ol\ = IQl]([€ l \) = \o L | T - Iff* ) M ( 2 ) 


where 

[ 0 /d kid = the mechanical stress vector 
(a/ ) 7 = the thermal stress vector 
{ < 7 / ) M == the moisture stress vector 

Next, Equation (2) is transformed from the material axes 1, 2, 3 to the laminated 
axes x, v and z (the z-axis coincides with the 3-axis) using the rotational transfor- 
mations as following: 


M = \T]-'\g,\ 
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where m = cos0, n = sin# and 6 = the angle between the I- and jc-axes. 
Thus, Equation (2) can be written in the form: 

M = l(?)|e| - Mr - Mm (3) 

Laminate Level Equations 

The strain-displacement relationship of the laminate at a distance z from the 
midplane of the laminate is: 


Id = Ml + Z(*i 


( 4 ) 


where 

|f| r = 

|f°| = is the strain vector at the midsurface of the laminate 

|&| r = \k,k y k, y | is the curvature vector at the midsurface of the laminate 

Substituting Equation (4) into Equation (3) we obtain the following: 

M = K?IU°I + ZlQ)\k\ - Mr - Mm ( 5 ) 

The vectors of the resultant forces and moments act at the midplane of the lami- 
nate with thickness h, and can be written in the following form: 


\N\ 

\m 



\a\dZ 


\o\zdz 


( 6 ) 

(7) 
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By substituting the ply stresses of the *th ply of Equation (5) into Equations (6) 
and (7) we finally obtain: 


N 


A 

B 


> 


N 


N 

A/J 


B 

D 


l*J 


[ A/J 

T 

\m[ 


where M J is the extensional stiffness matrix, [/?] is the coupling stiffness matrix 
and [D] is the bending stiffness matrix and are defined as following 


Ml = 

h/2 n 

\ Qldz = £ 

J h /2 K- \ 

e h* n 
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J***. ^K= 1 


where n is the number of the plies in the laminate. 

In order to perform the structural analysis of the fiber composite structure, it 
is required to know the mechanical and thermal properties of the laminate. 
Therefore the laminate is treated as homogeneous and possessing effective lami- 
nate engineering properties which are described in details in Reference [10]. 

Structural Level Equations 

The equations of equilibrium governing the linear dynamic response of the 
structure written in a matrix form are the following [8,9]: 

IM]\0\ 4- [D][U\ + [K]\U\ = [R(t)\ (9) 


where 

[A/J =the structural mass matrix 
[D] =the structural damping matrix 
[AT ] =the structural stiffness matrix 
(/?(/)( =the external applied load vector, as a function of time 
| V\ =the displacement vector 
| (7 ) =the velocity vector 
\U\ =the acceleration vector 


The equilibrium Equations (9), in the case of the free vibration problem are 
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simplified by setting the force vector equal to zero, \R(t)\ = 0. If the damping 
effect is not taking into account then the damping matrix is equal to zero, (D] = 
0. The equilibrium equations are simplified as following: 

[M}\0] + [K]\U\ = 0 (10) 


The structural matrices [M] and [K] are computed as follows: the continuous 
structure is discretized in a number of finite elements. If [H] is the shape function 
of a finite element of the structure then: 

(a) the mass matrix of an element is given by the equation: 


[m] 


q\ H\ T \H\ dV 


where q is the mass density of the element of the structure. 

(b) The stiffness matrix of an element is given by the equation: 

[*1 = J \B] T [C)[B\dV 

where [C] is the material matrix of the element, which contains the aniso- 
tropic effective properties of the laminate and is related to the stress-strain 
relationship as following: 


M = [C]| f | 

The matrix \B\ is related to the strain displacement relationship as follow- 
ing: 


|c| = lfl]l«l 

where |fl| = [L][N\, [ L\ is a linear operator and [WJ is the matrix which 
consists of the shape functions [//] [12]. 

The structural mass matrix [A/] is assembled from the element mass matrices 
[m]j for j = I, nel, where nel is the number of finite elements in the structure. 
Similarly, the structural stiffness matrix [Af] is assembled from the element 
stiffness matrices \k\j , j = 1, nel. 

After the computation of the [Af] and [AT] matrices the next step is to solve 
Equation (10). The solution of Equation (10) can be postulated to be of the follow- 
ing form: 


| (7 1 = Msin «(/ - r 0 ) 


(ID 
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Substituting Equation (11) into Equation (10) we obtain the following: 

[*1M = co 2 (12) 

For the computation of the vibration frequencies (eigenvalues) and the mode 
shapes (eigenvectors), CSTEM incorporates the following two methods: the de- 
terminant search and subspace iteration [9,11]. The search method is intended 
more for a small number of equations, while subspace iteration is more efficient 
for a large number of equations. 

GEOMETRY, LOAD HISTORY AND FINITE ELEMENT 
MODEL OF THE LAMINATED SHELLS 

The fiber composite thin shell had a cylindrical geometry with the following 
dimensions. The ratio of the longitudinal length L to the inner radius R was equal 
to L/R = 4. The ratio of the inner radius R to the thickness t of the shell was 
R/t — 33.3. The radius R was 25.4 cm. The boundary conditions were fixed at 
one end and free at the other (Figure 4). The laminate consists of continuous 
fibers made of the graphite material T300 embedded in intermediate modulus 
high-strength (IMHS) matrix. The material properties of the fiber and the matrix 
were taken from the data bank available in CSTEM (Tables 1 and 2, respec- 
tively). The fiber volume fraction (FVR) was 55 percent and the moisture content 
was 2 percent. A balanced and symmetric laminate was assumed with fiber orien- 
tations of[±0) 2 and plies of equal thickness. 

The simulation of the processing of the laminated cylinder was taken into ac- 
count. It started from the curing temperature T t of the matrix equal to 215.5 °C 
(Table 2) and gradually reduced until the oeprating temperature was equal to T g l2 
(Figure 4). At the end of the processing, residual stresses were induced in the 
constituent materials of the composite structure. 

For the free vibration analysis of the shell structure the effect of damping was 
neglected. The computational simulation was performed using three dimensional 


Table 1. T300 graphite fiber properties at room temperature 


Modulus in longitudinal direction. GPa 220.7 

Modulus in transverse direction, GPa 13 79 

In-plane Poisson s ratio 0.20 

Out-of-plane Poisson's ratio 0.25 

In-plane shear modulus, GPa 8 96 

Out-of-plane shear modulus, GPa 4.827 

Thermal expansion coefficients in longitudinal direction, 10 6 mm/(mm-°C) 0.99 

Thermal expansion coefficients in transverse direction, 106 mm/(mm-°C) 10.00 

Thermal conductivity in longitudinal direction, W/(m-K) 83.69 

Thermal conductivity in transverse direction, W/(m-K) 8.369 

Fiber tensile strength, GPa 2.413 

Fiber compressive s;rength, GPa 2.069 
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Table 2. IMHS matrix properties at room temperature 


Modulus in longitudinal direction, GPa 3.448 

Modulus in transverse direction, GPa 3.448 

In-plane Poisson’s ratio 0.35 

Out-of-plane Poisson’s ratio 0.35 

Thermal expansion coefficients in longitudinal direction, 10 6 mm/(mm-°C) 64.8 

Thermal expansion coefficients in transverse direction, 10 6 mm/(mm-°C) 64.8 

Thermal conductivity in longitudinal direction, W/(m-K) 2.16 

Thermal conductivity in transverse direction, W/(m-K) 2.16 

Matrix tensile strength, MPa 103.4 

Matrix compressive strength, MPa 241.3 

Matrix shear strength, MPa 89.6 

Void fraction 0 225 

Glass transition temperature, °C 215.55 


finite element analysis. One (eight node) element through the thickness of the 
composite shell structure was used. The finite element mesh consists of 760 
nodes and 360 elements (Figure 5). 

RESULTS AND DISCUSSION 

In this section results obtained for the different [±0] 3 angle-ply composite 
shells are presented and discussed. The parameters investigated include: 

1 . effects of fiber orientations 

2. effects of cylinder length 

3. effects of laminate thickness 

4. effects of temperature 

5. effects of fiber volume fraction 

6. effects of different ply thicknesses 

Effects of Cylinder Length 

The influence of the cylinder length in conjunction with the ply angle 6 was ex- 
amined to study the free vibration behavior of the composite shell. Results ob- 
tained for natural frequencies and mode shapes are summarized below. 

THE NATURAL FREQUENCIES 

Three thin composite shells with different length ratios L/R of 2, 4, and 6 were 
examined. The computed values of the first natural frequency is plotted versus the 
ply angle 6 and the three L/R ratios in Figure 6. The computed values are normal- 
ized with respect to the maximum frequency of 554.9 cps. Correspondingly, the 
values of the second frequency are normalized with respect to the maximum fre- 
quency of 609.6 cps and are plotted in Figure 7. Finally, the values of the third 
frequency are normalized with respect to the maximum frequency of 871.8 cps 
and are shown in Figure 8. The following conclusions are reached. 


Figure 5. Three dimensional finite element mesh and transverse cross sectional area of 
laminated shell structures. 


EFFECTS OF THE LENGTH RATIO L/R 
MAXIMUM FREQUENCY - 554.92 CPS 

-a- L/R-2 - t- L/R-4 - L/R-6 



Figure 6. Influence of the cylinder length on the lowest frequency of the laminated shell 
R/t = 33.3, FVR = 55%, moisture = 2%. 
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EFFECTS OF THE LENGTH RATIO l/R 
MAXIMUM FREQUENCY - 009.8 CPS 


-A- L/R-2 L/R-4 -a- l/R-0 



Figure 7. Influence of the cylinder length on the second frequency of the laminated shell. 


EFFECTS OF THE LENGTH RATIO L/R 
MAXIMUM FREQUENCY - 871.8 CPS 


“A- L/R-2 “ L/R-4 L/R-6 



PLY ANGLE, DEGREES 

Figure 8 . Influence of the cylinder length on the third frequency of the laminated shell. 
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Table 3. The percent difference in the natural frequencies between the 
short (L/R = 2) and the long (L/R = S) cylinders . 


Freq. 

to]. 

( ± 1 5|* 

I ±30]* 

[ ± 45]a 

[±60] 2 

t ± 75]j 

I ±90]* 

First 

57% 

60 7% 

63.5% 

70.7% 

73.3% 

74.1% 

75.4% 

Second 

60% 

59.6% 

67% 

65.9% 

63.5% 

55.7% 

40% 

Third 

58% 

45 4% 

50% 

55.4% 

54.9% 

57.7% 

57.3% 


Frequency ( L/R ~ 2) - Frequency ( UR = 6) 

Percent difference = * 100% 

Frequency (UR - 2) 


1. For all L/R , the natural frequencies increased for 0 < 0 < 22.5 and 
decreased for 22.5 < $ < 90. The maximum frequencies (first, second and 
third) occurred at 0 — 22.5. 

2. The shorter the shell the higher the frequencies as expected, since shorter 
shells are stiffer than the longer shells. 

3. The percent difference of the natural frequencies between the short 
(L/R = 2) and the long ( L/R = 6) shells for the different ply angles are tabu- 
lated in Table 3. It is seen that the percent difference of these natural frequen- 
cies is significant and dependent on the ply orientation angle. 

THE MODE SHAPES 

The behavior of the mode shapes of the composite shell with L/R = 6 and 

R/t = 33.3 for the first three frequencies are presented graphically as follows: 

I. First frequency (lowest frequency): the mode shape of the first frequency is 
shown in Figure 9. It is primarily bending indicating that the shell behaves 
like a cantilever beam. 



Figure 9. The first mode shape of the {±0 ] 2 angle-ply laminated shell. 
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Figure 10. The second mode shape of the [±0] 2 angle-ply laminated shell . 


2. Second frequency: the mode shape of the second frequency is shown in 
Figure 10. It is a combination of torsional and breathing (or radial) modes. 

3. Third frequency: the mode shape of the third frequency is shown in Figure 11. 
It is a combination of torsional, breathing (or radial) and axial (or longitudi- 
nal) modes. These mode shapes are summarized in Table 4 for convenience. 

Effects of the Laminate Thickness 

The influence of the laminate thickness of the thin composite shell versus the 
ply angle 0 for the three different laminate thickness ratios R/t of 20, 33.3 and 100 
are presented in Figure 12. The computed values are normalized with respect to 
the maximum frequency of 264.4 cps. The following conclusions were reached 
from Figure 12. 



Figure 11. The third mode shape of the [±$] 2 angle-ply laminated shell. 
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Table 4. Mode shapes of the laminated shell (LIR = 6, R/t = 33.3 
and fvr = 55%) as a function of the ply angle. 


Laminate 

First Frequency 

Second Frequency 

Third Frequency 

(0]. 

Beam mode 

Torsional and radial modes 

Torsional, axial and 
radial modes 

I±15J* 

Beam mode 

Torsional and radial modes 

Torsional, axial and 
radial modes 

(± 22 . 5)2 

Beam mode 

Torsional and radial modes 

Torsional, axial and 
radial modes 

[ ± 3012 

Beam mode 

Torsional and radial modes 

Torsional, axial and 
radial modes 

1 ± 45] 2 

Beam mode 

Torsional and radial modes 

Torsional, axial and 
radial modes 

l±60] 2 

Beam mode 

Torsional and radial modes 

Torsional, axial and 
radial modes 

[ ± 75)2 

Beam mode 

Torsional and radial modes 

Torsional, axial and 
radial modes 

[±90fc 

Beam mode 

Torsional and radial modes 

Torsional, axial and 
radial modes 


1. For all R/t, the first natural frequencies increased for 0 < 0 < 15 and 
decreased lor 15 < 0 < 90. The maximum frequency occurred at 0 = 15. 

2 . The thicker the composite shell the higher the values of the natural frequency. 
Thus thick composite shells (R/t = 20) had the highest frequencies, while 
thin composite shells (R/t = 100) had the lowest frequencies. 

3. The percent; difference of the natural frequencies between the thick laminated 
shell (R/t = 20) and the thin laminated shell (R/t = 100) are summarized in 
Table 5. As can be seen this difference is relative small. 

Therefore, the influence of the laminate thickness ratio R/t had a small effect 
on the first natural frequency of the laminated shell. 

Effects of Temperature 

The effects of the temperature profiles through the laminate thickness on the 
free vibration in conjunction with the ply angle 0 were examined. Three tempera- 
ture profiles with uniform moisture (2 percent) were considered: the first profile 
was uniform with temperature equal to room temperature (21 °C), the second 
profile was also uniform with temperature equal to TJ2 (107.7°C) and the third 
profile was varying linearly, with temperature equal to TJ2 in the inner surface 
and room temperature at the outer surface of the composite shell. The length 
ratio L/R was 4 and the laminate thickness ratio L/t was 33.3. The values of the 
first natural frequency are plotted versus ply angles 0 and temperature profiles in 
Figure 13. The computed values of the first natural frequency were normalized 
with respect to the maximum value of 276.8 cps. The following observation are 
made. 
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EFFECTS OF THE LAMINATE THICKNESS RATIO 
MAXIMUM FREQUENCY - 284.4 CPS 

R/t-20 R/t-33.3 -A- R/t-100 



Figure 12. influence of the laminate thickness ratio R/t on the lowest frequency of the lami- 
nated shell. UR * 4, FVR = 55% and moisture = 2%. 


1 . For all temperature profiles, the natural frequencies increased for 0 < $ < 
20 and decreased for 20 < 0 < 90. The frequency peaked at 0 = 20. 

2. The composite shell at room temperature had the highest frequencies, fol- 
lowed by the composite shell with the linear profile and finally the composite 
shell with the uniform temperature equal to TJ2. 

3. The percent difference in natural frequencies between the uniform profiles at 
room temperature and the T t !2 temperature is low (Table 6), and is almost in- 
dependent of the ply angle 0. 

Therefore, the temperature profiles had a small effect on the natural frequency of 

the composite thin shell throughout the angle-ply range. 


Table 5. The percent difference in the natural frequencies between the 
thick laminate (R/t = 20) and the thin laminate fR/t = 100). 


Freq. 

[01, 

[±15h 

l±30]. 

[±45), 

[±60], 

[±75], 

l±90h 

First 

2.35% 

2.5% 

4.5% 

9.62% 

1.8% 

0.8% 

1.2% 


Frequency {R/t = 20) - Frequency { R/t = 100) 

Percent difference = * 100% 

Frequency {RA = 20) 
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EFFECTS OF THE TEMPERATURE 
MAXIMUM FREQUENCY-276.8 CPS 




Figure 13. Influence of the temperature profile on the lowest frequency of the laminated 
shell. UR = 4, Fl/t = 33.3, FVR = 55%, moisture *= 2%. 


Effects of Fiber Volume Fraction 

The effects of the fiber volume fraction in conjunction with the ply angle 6 on 
the free vibration behavior of the composite shell were evaluated. Three fiber 
volume fractions were examined: 55, 60 and 65 percent. The length ratio L/R was 
4 and the laminate thickness ratio R/t was 33.3. The values of the first natural 
frequency are plotted for the different ply angles 8 and for the various fiber 
volume fractions in Figure 14. The values of the natural frequency are normalized 
with respect the maximum value of 278.6 cps. The following conclusions were 
reached. 


Table 6. The percent difference in the natural frequencies between the 
laminate shelis in room temperature ( T = 21 C) and 
In high temperature (T = Tq/2). 


Freq. 

[0I« 

[±15], 

[ ± 30] 2 

[±45], 

(±60], 

[±751, 

l±90J, 

First 

3.8% 

2.5% 

3.1% 

3.6% 

3.3% 

3% 

3.5% 


Frequency If * 21 C) - Frequency (7 * TqI 2) 


Percent difference = 


Frequency (T ■ 21 C) 


* 100 % 
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EFFECTS OF THE FIBER VOLUME FRACTION 
MAXIMUM FREQUENCY-270.0 CPS 

FVR-55% FVR-00% FVR-05% 



Figure 14. Influence of the fiber volume fraction on the lowest frequency of the laminated 
shell. L/R = 4, R/t =* 33.3, moisture = 2%. 


1 For all fiber volume fractions the natural frequencies increased for 0 < 0 < 
20 and decreased for 20 < $ < 90. The maximum natural frequency occur- 
red at 0 - 20, which is the same for all the cases investigated. 

2. The composite shell with the fiber volume fraction of 65 percent had higher 
frequencies than the shell with the fiber volume fraction of 60 percent which 
had higher frequencies than the shell with fiber volume fraction of 55. 

3. The percent difference in the natural frequencies between the two laminates 
with fiber volume fraction of 65 and 55 is very low and is shown in T&ble 7. 

Therefore, the influence of the fiber volume fraction had small effect on the 

natural frequency of the composite shell. 


Table 7. The percent difference in the natural frequencies between the 
laminates with fvr of 65°A> and fvr of 55%. 


Freq. 

[0]4 

[±15Jj 

f±30] a 

l±45h 

[±60fe 

1*751, 

[±90J, 

First 

6.06% 

6.05% 

6.1% 

6.3% 

6.2% 

5.4% 

5.8% 


Percent difference = 


Frequency (fvr = 65%) - Frequency (fvr 5 65%) 
* 100 % 


Frequency {RA f = 20) 
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EFFECTS OF THE PLY THICKNESS 
NORM MIZED MAXIMUM FREQUENCY -262.3 CPS 



Tie] 

U*l 2 

__ 0 __ (3t/8.t/8.t/8.3t/8) 
[i0] 2 

— V — { t/4. t/4,t/4. t/4) 

t±ei 4 

D ... (t/8,t/8 t/8) 


Figure 15. Influence of the ply thickness on the lowest frequency of the laminated shell. 
L/R - 4, R/t - 33.3, FVR - 55%, and moisture = 2%. 


Effects of Ply Thickness 

The effects of laminates with various ply thicknesses, in conjunction with the 
ply angle 6 , were examined on natural frequencies of the composite shell. The 
length ratio L/R was 4, the laminate thickness ratio R/t was 33.3 and the fiber 
volume fraction was 55 percent. Four laminated shells with different ply thick- 
nesses and the same laminate thickness (f) were examined: (1) a [ ±0] angle-ply 
laminate with equal ply thicknesses (tf 2, t! 2); (2) a [±0] 2 angle-ply laminate 
with nonuniform ply thicknesses (3r/8, tf 8, tf 8, 3f/8); (3) a [dt0] 2 angle-ply lami- 
nate with uniform ply thicknesses (f/4, // 4, f/4, f/4); (4) a [ ±d] 4 angle-ply lami- 
nate with equal ply thicknesses (f/8, tf 8, ...» t/8). The computed values of the 
first natural frequency of these examined laminates are plotted for different ply 
angle 0 in Figure 15. The following conclusions were reached. 

1 . The natural frequencies increased for 0 < 0 < 20, decreased for 20 < 6 < 
90 and peaked at 0 = 20. 

2. The laminate with different ply thicknesses had insignificant effect on the 
natural freouencies of the composite thin shell. 

SUMMARY OF THE RESULTS 

A computational simulation of the free vibration of hot, wet composite, thin 
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shells was performed by using a three-dimensional finite element analysis com- 
puter code. 

The simulation included the effects of parameters such as ply angle, cylinder 
length, laminate thickness, temperature, fiber volume fraction and ply thickness 
on the free vibration behavior of the composite shells. The important results are 
summarized as follows. 

1 The fiber orientations and the length of the laminated shell significantly in- 
fluenced the free vibration behavior (frequencies and mode shapes) of the 
composite shell. 

2. The fiber volume fraction, the laminate thickness and the temperature profile 
through the laminate thickness had a small effect on the free vibration be- 
havior of the composite shell. 

3. The laminates with different ply thicknesses had an insignificant effect on the 
free vibration behavior of the composite shell. 

4. A single, through -the-thickness, eight-node, three-dimensional, composite 
finite element analysis appears to be sufficient for investigating the free vibra- 
tion behavior of thin composite angle-ply cylindrical shells. 
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